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Abstract. In this article, we present a mathematical theory of the Walker 
circulation of the large-scale atmosphere over the tropics. This study leads to 
a new metastable state oscillation theory for the El Nino Southern Oscillation 
(ENSO), a typical inter-annual climate low frequency oscillation. The mathe- 
matical analysis is based on 1) the dynamic transition theory, 2) the geometric 
theory of incompressible flows, and 3) the scaling law for proper effect of the 
turbulent friction terms, developed recently by the authors. 



1. Introduction 

The atmosphere and ocean around the earth are rotating geophysical fluids, 
which are also two important components of the climate system. The phenomena 
of the atmosphere and ocean are extremely rich in their organization and complex- 
ity, involving a broad range of temporal and spatial scales. According to J. von 
Neumann |12j . the motion of the atmosphere can be divided into three categories 
depending on the time scale of the prediction. They are motions corresponding 
respectively to the short time, medium range and long term behavior of the at- 
mosphere. One of the primary goals in climate dynamics is to document, through 
careful theoretical and numerical studies, the presence of climate low frequency 
variability, to verify the robustness of this variability's characteristics to changes 
in model parameters, and to help explain its physical mechanisms. The thorough 
understanding of this variability is a challenging problem with important practical 
implications for geophysical efforts to quantify predictability, analyze error growth 
in dynamical models, and develop efficient forecast methods. One typical source of 
the inter-annual climate low frequency variability is El Nino-Southern Oscillation 
(ENSO) , which is associated with the Walker circulation of the atmosphere over the 
tropics. The main objective of this article is to introduce a mathematical theory of 
the general circulation of the atmosphere over the tropics. 

The modeling part of the study is based on two important ingredients. First, 
from the general circulation point of view, there are three invariant regions of the 
atmosphere: the northern hemisphere, the southern hemisphere and the region 
over the tropics corresponding to zero latitude. The second is the new scaling law 
introduced recently by the authors, leading to correct circulation length scales. 

The model is then analyzed using the dynamic transition theory and the geo- 
metric theory of incompressible flows, both developed recently by the authors. We 
refer the interested readers to |H] and the references therein for the geometric theory, 
and to the appendix, [4j [10] and the references therein for the dynamic transition 
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theory. We mention in particular here that the main philosophy of the dynamic 
transition theory is to search for the full set of transition states, giving a complete 
characterization on stability and transition. This complete set of transition states 
represent the "physical reality" after the transition, and are described by a lo- 
cal attractor, rather than some steady states or periodic solutions or other type of 
orbits as part of this local attractor. With this theory, many longstanding phase 
transition problems are either solved or become more accessible leading to a num- 
ber of physical predictions. For example, the study of phase transitions of helium-3 
leads not only to a theoretical understanding of the phase transitions to superfluid- 
ity observed by experiments, but also to such physical predictions as the existence 
of a new superfluid phase C for liquid helium-3 [9J . 

The main results obtained in this article, on the one hand, verifies the general 
circulation patterns over the tropics associated with the Walker circulation, and, on 
the other hand, lead to a new mechanism of El Nino Southern Oscillation (ENSO), 
which will be published in an accompanying paper |5]. This new mechanism of 
the ENSO amounts to saying that ENSO, as a self-organizing and self-excitation 
system, with two highly coupled processes. The first is the oscillation between 
the two metastable warm (El Nino phase) and cold events (La Nina phase), and 
the second is the spatiotemporal oscillation of the sea surface temperature (SST) 
field. The interplay between these two processes gives rises the climate variability 
associated with the ENSO, leads to both the random and deterministic features 
of the ENSO, and defines a new natural feedback mechanism, which drives the 
sporadic oscillation of the ENSO. 

The paper is organized as follows. In Section 2, we introduce the atmospheric 
circulation model over the tropics, which is analyzed in Sections 3 and 4, with 
concluding remarks given in Section 5. 



Physical laws governing the motion and states of the atmosphere and ocean can 
be described by the general equations of hydrodynamics and thermodynamics. As 
discussed in the Introduction, the atmospheric motion equations over the tropics 
are the Boussinesq equations restricted on 9 — 0, where the meridional velocity 
component uq is set to zero, and the effect of the turbulent friction is taking into 
considering using the scaling law derived in [ID] : 



2. Atmospheric Model over the Tropics 
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Here Oi — Cih 2 {i — 0,1) represent the turbulent friction, a is the radius of the 
earth, the space domain is taken as M = S\ x (a, a + h) with S\ being the one- 
dimensional circle with radius a, and 



u v d d 
a 0(p oz 



1 d 2 



dz 2 



a 2 dtp 2 
For simplicity, we denote 

(xi,x 2 ) = (cup, z), (ui,u 2 ) = (u v ,u z ). 

In atmospheric physics, the temperature T\ at the tropopause z = a + h is a 
constant. We take To as the average on the lower surface z = a. To make the 
nondimensional form, let 

x — hx' , t — h 2 t' / k, u — kv! /h, 

T = (To - Ti)T' + To - (T - T )x' 2 , 

P = Kvpop'/h 2 - gp (hx' 2 + a T (T - T 1 )hx' 2 /2). 

Also, we define the Rayleigh number, the Prandtl number and the scaling laws by 

UTg(T - Tx)h 3 „ v 



(2.2) 



R 



Pr 



Omitting the primes, the nondimensional form of (2.1 1 reads 
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dT 
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where (xi,x 2 ) £ M = (0,27rr ) x (r ,r + 1), Sq and 8± are as in (2.2), (u • V) and 
A as usual differential operators, and 

(2.4) uj= . 

K 

The problem is supplemented with the natural periodic boundary condition in 
the Xi-direction, and the free-slip boundary condition on the top and bottom bound- 
ary: 



(2.5) 
(2.6) 



(u,T)(x 1 +2nr ,x 2 ,t) = (u, T)(x 1 , x 2 , t) 
dui 



u 2 = 0, 



dx 2 



0,T = ip(xi) 



u 2 = 0, = 0, T = 
ax 2 



at x 2 = r , 
at x 2 = r + 1 . 
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Here <p(x\) is the temperature deviation from the average To on the equatorial 
surface and is periodic, i.e., 

2-Kro 

ip{x\)dxi — and <p{x\) = <p(x\ + 27rr ). 

The deviation tp(xi) is mainly caused by a difference in the specific heat capacities 
between the sea water and land. 



3. Walker Circulation under the Idealized Conditions 

In an idealized case, the temperature deviation ip vanishes. In this case, the 
study of transition of (2.3 1 is of special importance to understand the longitudinal 
circulation. Here, we are devoted to discuss the dynamic bifurcation of (2.3), the 
Walker cell structure of bifurcated solutions, and the convection scale under the 
idealized boundary condition 

tp(x\) = for any < x\ < 2irro, 
Q = in pT3l. 



(3.1) 



For the problem ( |2.3| ) with ( |2.5[ ) and ( |2.6| , let 

H = {(u,T) e L 2 (Mf | divu = 0, (u,T) satisfies fiDfy, 



Hi = {(u,T) e H 2 {Mf n H | (u,T) satisfies one of (j^) and (|2j|)}. 
Then, define the operators Lr = A + Br and G : Hi — > H by 



Ayj = P(P k (A Ul 
(3.2) B R ijj = P{-uu 2 , w«i + Pr RT, u 2 ), 



— ~ o M i): Pr ( Au 2 7, 0^2), AT), 

ro ox i ro ox± 



= P(((« • V)«i + — ), ((« • V) W2 - (u • V)T), 

where ip = (u,T) E Hi, P : L 2 (M) 3 -> ff is the Leray Projection, and 5 ■ = ^ + $ 

(i = 0, 1). Under the definitions (3.2), the problem (2.3 1 with one of (2.5 1 and (2.6) 
is equivalent to the following abstract equation 

(3.3) ^W + GW- 
Consider the eigenvalue problem 

(3.4) L fl V = P(K)if>, 
which is equivalent to 

dp 
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(3.5) 



Au 2 
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at + u 2 — pr. 

divu = 0, 



- S[u 2 
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dp 



RT--^ = P-'(3u 2 , 
axo 



with the boundary conditions (2.51 and (2.6 1. 
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We shall see later that these equations (3.5) are symmetric, which implies, in 
particular, that all eigenvalues (3j(R) are real, and there exists a number R c , called 
the first critical Rayleigh number, such that 



(3.6) 



(3.7) 



r < if R < R c , 

Pi(R)l -0 if R = R C , 

{ > if R > R c 

Pj{R c ) <0 



for i = 1,2, 



for j > 2. 



The following theorem provides a theoretical basis to understand the equatorial 
Walker circulation. 



Theorem 3.1. Under the idealized condition (3.1), the problem (2.3) with (2.5) 



(2) 



and (2.6) undergoes a Type-I transition at the critical Rayleigh number R = R c . 
More precisely, the following statements hold true: 

(1) When the Rayleigh number R < R c , the equilibrium solution (u, T) = is 
globally stable in the phase space H. 

When R c < R < R c + 5 for some S > 0, this problem bifurcates from 
((it, T), R) = (0,i? c ) to an attractor Ar = S , consisting of steady state 
solutions, which attracts H\T , where T is the stable manifold of '(u,T) = 
with codimension two. 

(3) For each steady state solution ipR — (ur,Tr) £ Ar,Ur is topologically 
equivalent to the structure as shown in Figure \3.1\ 

(4) For any initial value tpo = (mo,?o) £ H \ (T U E), there exists a time 
tg > such that for any t > to the velocity field u(t, tpo) is topologically 
equivalent to the structure as shown in Figure 3.2 either (a) or (b), where 
tp = (u(t,ipo),T(t,ipo)) is the solutions of the problem with ip(0) — ipo, and 

ro + l 



E = {(u,T)eH\ / Ul dx 2 = 0} 




r„+ 1 



Figure 3.1. The cell structure of the steady state solutions in the 
bifurcated attractor Ar. 



Proof of Theorem\3.1\ We proceed in the following several steps. 



Step 1. We shows that equations (2.3 1 have an equivalent form as the classical 
2D Benard problem. 

Since the velocity field u defined on M = S 1 x (ro,ro + 1) is divergence-free, 
there exists a stream function tp such that 

dtp dp 

OX2 OX 1 
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(a) (b) 



Figure 3.2. The Walker cell structure with the cells separated by 
a cross channel flow: (a) a westbound cross channel flow, and (b) 
an eastbound cross channel flow. This cross-channel flow pattern 
has the same topological structure as the Walker circulation over 
the tropics and the Branstator-Kushnir waves in the atmospheric 
dynamics [2I2]- 



satisfying the given boundary conditions. Therefore, the following two vector fields 

(-W!l2,W!il) = wVjJ, 



(3.8) 



^ 2 du 2 2 du\ 2 ^ dip 

vq dxi ' ro dx\ rp dx\ 



arc gradient fields, which can be balanced by Vp in ( 2.3 ) . Hence, ( 2.3 1 are equivalent 
to 



— + (u ■ VjMi H = Pr 

at r Q 



(3.9) 



du 2 

~dt 

dT 

St ' 
divu 



+ (u ■ V)u 2 = Pr 



(u ■ V)T = AT - 
= 0. 



Aui — SqUi 
Au 2 - S[u 2 + RT 



dp 
dxi 



dp 
dx 2 



u 2 , 



Therefore, equation (3.3 1 is also an abstract form of (3.9 1. Thus, equations 



(3.4 1 and (3.5) are symmetric. It is clear that the operator G defined in (3.2) 
is orthogonal. Hence, the attractor bifurcation theorem for the Rayleigh-Benard 
convection proved in [5] is also valid for the problem (2.3 1 with (2.5 1, ( |2.6[ ). Thus, 
this transition of the problem is Type-I (i.e., continuous), and Assertion (1) is 
proved. 

Step 2. Proof of Assertion (2). We only need to verify that the attractor 
Ar is homeomorphic to a circle, and consisting of singular points of (3.3 1. We 
consider the first eigenvectors of ( 3 . 5 1 at R = R c . By (3.8|, equations (3.51 at 
R c are equivalent to the form 



R 



(3.10) 



Aui — S' ui 



dp 
dx\ 

Au 2 - 5[u 2 + R C T 

AT + u 2 = 0, 
divu = 0, 



0. 

dp 
dx 2 



0. 



with the boundary conditions (2.51 and (2.6 1 
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As in j3j Ej, we see that the multiplicity of R c in (3.101 is m = 2, and the 
corresponding eigenvectors are given by 

... . ( kx\ i , . k . l kx\ , l \ 
(3.11) = {ui,u 2 ,T) = - sin ii (x 2 ), — cos iitehcos <P(:E2) , 



(3.12) Vi 



ui,u 2 , T) = (cos — — -H'(x2), — sin — -H(x2),ain — — -$(x2) ) , 
ro r r r 



(3.13) 



where the functions H(x 2 ) and Q(x 2 ) satisfy that 

' (L> 4 - (a 2 + a 2 ,)!) 2 + a 2 a\)H = aR c <S>, 
- (D 2 - a 2 )$ = aH, 

(3.14) H = H" = 0, $ = at a; 2 = r , r + 1, 

where D = d/dx 2 ,a = k/r , a 2 — a 2 + 8[ (i = 0, 1). 

Since (3.131 with ( |3.14[ ) are symmetric, we define a number a by 

(3.15) a= _?_<G(* > ^) + G(^i ) *),V'i>, 



101 



where a is the number as in (A.6), and \& is the center manifold function and by 
[4] it satisfies that 

(3.16) — LrJ5> = G(ipx,ipi) + higher order terms, 

and Lr is defined by (3.2 1, and G(tj},ip) defined by 

G$, 4>) = P ((u • V)«i + — , (u • V)u 2 + — , (u • V)T 



for any ip = (u,T),ip = («, T) € Obviously, for any 0, G ffi, we have 
<G(^),^>= 0, 
<G(0,0),0>=-<G(0>),0>- 



Hence, we derive from (3.15) and (3.161 that 
' <G(Vi,Vi).*>=- 



1 



ll^il 



IV'ill 2 



Due to (3.8 1, — is a symmetric sectorial operator. Thus, we have 

' Ki-Lny/^^-L^) 1 ^ ><0. 



IV'ill 2 



Then, as in [6], Assertion (2) follows. 

Step 3. Proof of Assertion (3). By Step 1 and the attractor bifurcation 
theorem, each steady state solution <p S Ar can be expressed as 

(3.17) V = |/3 1 (i?)/a| 1 / 2 (^ 1 + #x) + o(Jh(R)), 

where ipx,ipi are given by (3.111 and (3.121, x 2 + y 2 = 1, and Pi(R) is the first 
eigenvalue of (3.4 1 satisfying (3.6 1. 

Let if = (e, T). Then, by (3.17) the velocity field e is written as 

r(R) cos — (si + e)H'{x 2 ) + o(r(i?)), 
— rfifc) sin — (n + 9)H(x 2 ) + oMR)), 
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where r(R) = \(3i (R)/a \ 1/2 , (x, y) = (cos 0, sin 6). 

By (3.131 and (3.14), H(x 2 ) = sinjir(x 2 — r ),j G Z. But, for the first eigenvec- 
tors ip\ and tpi, we have 



Thus, e is expressed as 
(3.18) 



H(x 2 ) = sin7r(x 2 - r ). 



e + o(r(i?)), 



(3.19) 



eo 



irr(R) cos — (xi + 0) cos7r(x2 — r ), 
— r(i?) sin — (a?i + 6) sin 7t(;e2 — r ), 



It is easy to check that eo is regular, therefore e is also regular as R c < R < R c + 5 
for some S > 0. Obviously, the vector field eo is topologically equivalent to the 
structure as shown in Figure |3.1| 



By the structural stability theorem in [5], the vector field eo in (3.19) is not 
structurally stable in Hi, because the boundary saddle points on x 2 = Tq are con- 
nected to saddle points on x% = ro + 1, a different connected component. However, 
e is structurally stable in the space 

/•2irr /To + 1 



H = {(u,T) e H\ 



udx = 0}. 



To see this, we know that if u G H, then u = {ui, U2) ^constant, and Ui has the 
Fourier expansion 



It follows that 
(3.20) 



ui = y^(xfcj ■ 

r +l 



\-y kj sin )cosj7r(a;2 -r J. 



uidx 2 = Vu = "2) G H 



Accor ding to the Connection Lemma (Lemma 2.3.1 in Ma and Wang [S]), from 
(3.201 one can infer that a vector field u £ H is structurally stable if and only if 
u is regular, and all interior saddle points of tt are self-connected. Therefore, eo is 
structurally stable in H . 

If we can prove that the vector field e given by (3.181 is in H, then, as R c < R < 
R c + 5, e is topologically equivalent to eo. Hence, to prove Assertion (3), it suffices 
to verify that e E H. 

Obviously, H is an invariant space for the operator Lr + G defined by (3.2), and 
the orthogonal complementary H 1 - of H in H is 

H 1 - = {{u, T) £ H | u = (c, 0), T = 0, c G R 1 } 



It is readily to prove that all steady state solutions of ( 3.3 1 are in H. Thus, Assertion 
(3) is proved. 

Step 4. Proof of Assertion (4). For any initial value Vo = (uo,T ) G 
H\H,ipQ can be written as 

^o = cq + $o,$o G H, 

c = {a 0> 0,0)€H ± , a ^0. 
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Since G{ijj) G H^tp S H, the solution ip(t,ipo) of (3.3 1 with tp(0, i/j ) = ipo must be 
in the form 

In addition, we know that for each -0o G i?/r there is a steady state solution 
ip = (e, T) G yl_R such that 



lim 



(i, V'o) - vile = for any r > 1. 



Then, for any ip £ H \ (H U T) , it follows that there exists a time > such that 
the velocity field u(t,ipo) of tp(to,ipo) is topologically equivalent to the following 
vector field for any t > to 



u — e + Coe 



e+(a o e- ,5 « t ,0,0) 



where e is as in (3.18). 

The vector field e has the roll structure as shown in Figure [3~lj and c e~ s o t is a 
parallel channel flow. Hence, it is easy to show that u — e + coe~ Sot is topologically 
equivalent to the structure as shown in Figure [3~2] (a) as ao < 0, and to the structure 
as shown in Figure 3.2 (b) as ao > 0, for t > t large. Assertion (4) is proved. □ 



4. Walker circulation under natural conditions 
We now return the natural boundary condition 

ip(xx) ^ and Q^O. 



In this case, equations (2.3) admits a steady state solution 

(4.1) i> = {V,J), (V=(Vi,Va)). 

Consider the deviation from this basic state: 

u^u + V, T -^T+J. 



Then (2.3 1 becomes 
dui 



dt r a 
- (u • V)tn - (it • V)Vi 



Pr [Ami - S' U! - 

Vi V 2 
— u 2 til, 



dp 1 



(4.2) 



^ + (u ■ V) M2 - ^ = Pr [Au 2 - S[u 2 + RT-^L] 
at r ox 2 



dT 



2V\ 

(it ■ V)u 2 - (u ■ V)V2 + — ui, 
- (it • V)T = AT+ U2 - (it • V) J - (it • V)T, 



at 

di vit = 0. 

The boundary conditions are the free-free boundary conditions given by 
(it, T) is periodic in x\ — direction, 



(4.3) 



T = 0, u 2 = 0, - — = at x 2 = ro, r + 1. 
dx 2 
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Let the operators Lr = A+Br and G : Hi 
be denned by 



P(-(V • V)wi - (w • V)Vi 



-> H be as in (3.2 1, and C £ : Hi 
Vi 



H 



ro 



-u 2 



r 



-Ml, 



2Vi 



-(V ■ V)u 2 - (u ■ V)V 2 + —ui,-(u • V) J - 



(V ■ V)T). 



Then, the problem (4.2) and (|4.3| is equivalent to the abstract form 
(4.4) 



dv 



dt 



= Lr^ + £ £ iP + G(iP). 



Consider the eigenvalue problem 

(4.5) LriP + £^ = f3 e (R)^. 

It is known that \tp(xi)\ and Q ^ are small, with AT = T - T x = 100°C as 
unit. Hence, the steady state solution (V, J) is also small: 

\\(V,J)\\ L2 =e<^l. 



Thus, ( |4.2| is a perturbation equation of (2.3 1. 

Since perturbation terms involving ( V, J) are not invariant under the zonal trans- 
lation (in the xi-direction), for general small functions <p(xi) ^ 0, the first eigen- 
values of (4.2) are (real or complex) simple, and by the perturbation theorems in 
[4], all eigenvalues of linearized equation of (4.2 1 satisfy the following principle of 
exchange of stability (PES): 



for any 1 < i < m, 



for R < R e c , 
for R = R S C , 
for R > R e c , 

for any j > m + 1 , 

where m = 1 as (3f(R) is real, m = 2 as f3f(R) is complex near i?^, and i?J is the 
critical Rayleigh number of perturbed system (4.2 1. 




The following two theorems follow directly from Theorem 3.1 and the perturba- 
tion theorems in the appendix or in [4]. 



Theorem 4.1. Let Pi(R) near R = R 6 C be a real eigenvalue. Then the system 
has a transition at R = R £ c , which is either mixed (Type-Ill) or continuous 
(Type-I), depending on the temperature deviation <p(xi) 



Moreover, we have the 
following assertions: 

(1) // the transition is Type-I, then as R e c < R < R E C + 5 for some S > 0, the 
system bifurcates at R e c to exactly two steady state solutions ipi and ip 2 in 
H , which are attractors. In particular, space H can be decomposed into two 
open sets Ui,U 2 : 



H=Ui 



U 2 



UiHU 2 = 



ifj = 0EdUir\dU 2 , 



such that ipi € ?7, (i = 1, 2), and tpi attracts Ui. 
(2) // the transition is Type-Ill, then there is a saddle-node bifurcation at R = 
R* with R* < R E C such that the following statements hold true: 
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(a) if R* < R < R e c + 5 with R =/= R e c , the system has two steady state 
solutions if/ft and ip^ which are attractors, as shown in Figure J^.l 
such that 

= for R* <R< R%. 

(b) There is an open set U C H with G U which can be decomposed into 
two disjoint open sets U — UP + U 1 } with ip^ £ U± and ijr^ attracts 

(3) For any initial value ipo — ( u o, To) € for R> R* or ipo — (uq, Tq) E U+ 
for R > R e c , there exists a time t^>Q such that for any t > to the velocity 
field u(t,ipo) is topologically equivalent to the structure as shown in Figure 



3.2 either (a) or (b), where ip = (u(t,^o),T(t,^o)) is the solutions of the 



problem with ip(0) — ipo- 




ty R El Nino phase 



R Rayleign number 
ip B " La Nina phase 

Normal phase 

Figure 4.1. The transition diagram for the equatorial motion 
equations (2.3 1, R* is the saddle-node bifurcation point, R e c is the 
first critical Rayleigh number. 



Theorem 4.2. Let /3f(R) be complex near R = R e c . Then the system (4-2) bifur- 
cates from (ip,R) — (0,R E C ) to a periodic solution ipR(t) on R e c < R, which is an 
attractor, and ipR(t) can be expressed as 

(4.6) i> R (t) = A fl (cosp#i + sinpt^i) + o(\A R \,e), 

where Ar = a{R — R £ c ), a and p are constants depending on (p(xi), and tpi, ipi are 



first eigenf unctions of linearized equations of (2.3) 



5. Concluding Remarks 

In this article, a careful examination of the dynamic transitions and stability 
of the large-scale atmospheric flows over the tropics, associated with the Walker 
circulation and the ENSO, is given. The analysis and the results obtained show the 
following from the physical point of view. 



First, Theorems 3.1 4.1 and 4.2 provide the possible dynamical behaviors for 

Theorem 13.11 is for the ideal case 



the atmospheric circulation over the tropics, 
where the surface temperature profile is given as a constant, leading to a translation 
oscillation. Apparently, this result does not represent a realistic explanation to the 
ENSO. 
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Second, the time-periodic oscillation obtained in Theorem |4.2| does not repre- 
sent the typical oscillation in the ENSO phenomena either, as the Walker circulation 
does not obey the zonal translational oscillation given by the periodic solutions. 



Third, Theorem 4.1 does give a characterization of the ENSO, leading a correct 
oscillation mechanism of the ENSO between two metastable El Nino and La Nina 
events, which is further studied in 18 . 



Appendix A. Dynamic Transition Theory for Nonlinear Systems 

In this appendix we recall some basic elements of the dynamic transition theory 
developed by the authors |H[7], which are used to carry out the dynamic transition 
analysis for the binary systems in this article. 

A.l. New classification scheme. Let X and X\ be two Banach spaces, and 
X\ C X a compact and dense inclusion. In this chapter, we always consider the 
following nonlinear evolution equations 

(A.l) Tt =L x u + G(u,\), 

u(0) = (p, 

where u : [0, oo) — > X is unknown function, and A £ K 1 is the system parameter. 

Assume that L\ : X\ — > X is a parameterized linear completely continuous field 
depending continuously on A £ I 1 , which satisfies 

L\ = —A + B\ a sectorial operator, 
(A. 2) A : X\ — > X a linear homeomorphism, 

B\ : X\ — > X a linear compact operator. 

In this case, we can define the fractional order spaces X a for a € R 1 . Then we also 
assume that £?(•, A) : X a — > X is C r (r > 1) bounded mapping for some < a < 1, 
depending continuously on A € K 1 , and 

(A.3) G(u,X) = o(\\u\\ Xa ), VA e R 1 . 



Hereafter we always assume the conditions ( |A.2 1 and (A.3 1, which represent that 



the system (A.l I has a dissipative structure. 



Let the eigenvalues (counting multiplicity) of L\ be given by 

{&(A)eC | J = 1, 2, • • - } 

Assume that 

(A.4) Re /3 4 (A) { = if A = Ao, VI < % < m, 




(A.5) Re P 3 (Xq) < Vj > m + 1. 



The following theorem is a basic principle of transitions from equilibrium states, 
which provides sufficient conditions and a basic classification for transitions of non- 
linear dissipative systems. This theorem is a direct consequence of the center man- 
ifold theorems and the stable manifold theorems; we omit the proof. 
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Theorem A.l. Let the conditions (A.4\ ) and (A.5\ ) hold true. Then, the system 



(A.l) must have a transition from (u, A) = (0, Ao), and there is a neighborhood 
U C X of u = such that the transition is one of the following three types: 

(1) CONTINUOUS TRANSITION: there exists an open and dense set U\ C U 



such that for any <p 6 U\, the solution U\(t, p) of (A.l) satis fi 



k s 



lim limsup ||ua(£, <p)\\x = 0. 
A — >Ao t — >oo 



In particular, the attractor bifurcation of (A.l) at (0, Ao) is a continuous 
transition. 

(2) Jump TRANSITION: for any A < A < A + e with some e > 0, there is an 
open and dense set U\ C U such that for any if € U\, 

limsup \\u\(t, ip)\\x > 5 > 0, 

t — >oo 

where S > is independent of A. This type of transition is also called the 
discontinuous transition. 

(3) Mixed Transition: for any A < A < A + e with some e > 0, U can be 
decomposed into two open sets U± and (Uf not necessarily connected): 

U=U? + C/ 2 \ u£ n t/ 2 A = 0, 

such that 



lim limsup | u(t , ip) \\ x = V<y9 € U{ 

A^Ao t-^oc 

limsup \\u(t, ip)\\x > S > y<p e U. 

t — >oo 



A 



With this theorem in our disposal, we are in position to give a new dynamic 
classification scheme for dynamic phase transitions. 

Definition A.l (Dynamic Classification of Phase Transition). The phase transi- 



tions for (A.l) at A = Ao is classified using their dynamic properties: continuous, 
jump, and mixed as given in Theorem \A.1\ which are called Type-I, Type-II and 
Type-Ill respectively. 

An important aspect of the transition theory is to determine which of the three 
types of transitions given by Theorem | A . 1 1 occurs in a specific problem. Hereafter 
we present a few theorems in this theory to be used in this article, and we refer 
interested readers to |10j for a complete description of the theory. 



A. 2. Transitions from simple eigenvalues. We consider the transition of ( A.l ) 
from a simple critical eigenvalue. Let the eigenvalues (3j(X) of L\ satisfy (A.4) and 



( A.5[ ) with m = 1. Then the first eigenvalue /3i(A) must be a real eigenvalue. Let 



ei(A) and e*(A) be the eigenvectors of L\ and L* x respectively corresponding to 
/3i(A) with 

L Ao ei=0, L* x el = 0, <e 1 ,e* 1 >=l. 



Let Q(x, A) be the center manifold function of (A.l I near A = Ao- We assume that 
(A.6) < G(xex + $(sg, A ), A ), e* >= ax k + o(\x\ k ), 

where k > 2 an integer and a/Oa real number. 



14 



MA AND WANG 













"l 1 


u = 


V 2 



(a) 



(b) 



Figure A.l. Topological structure of the jump transition of (A.ll 
when A:=odd and a > 0: (a) A < Ao; (b) A > A . Here the 
horizontal line represents the center manifold. 



u = 



(a) 



(b) 



Figure A. 2. Topological structure of the continuous transition of 



(A.l I when A:=odd and a < 0: (a) A < A ; (b) A > \ . 



(a) 



(b) 



(c) 



Figure A. 3. Topological structure of the mixing transition of 

(a) A < A ; (b) A 



( |A.1| when fc=even and a ^ 0: (a) A < Ao; (b) A = Ao; (c) 
A > Ao- Here [/A is the unstable domain, and C/^ the stable do- 
main. 



Theorem A. 2. Assume (A. 4) and (A. 5) with m = 1, and (A. 6). If k=odd and 
a =/= in (A. 6) then the following assertions hold true: 

(1) If a > 0, then (A.l) has a jump transition from (0, Ao), and bifurcates on 
A < Ao to exactly two saddle points Vi and v% with the Morse index one, 

H3 



as shown in Figure 



(2) lfa<0, then (A.l) has a continuous transition from (0,Aq), which is an 



attractor bifurcation as shown in Figure A. 2 
(3) The bifurcated singular points and v% in the above cases can be expressed 
in the following form 



l/fe-i 



ei 



(A) + (|/3 1 | 1 / fc - 1 ). 
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Theorem A. 3. Assume (A. 4) and (A. 5) with m = 1, and (A. 6). If k= even and 
a =/= 0, then we have the following assertions: 

(1) (A.l) has a mixed transition from (0,Ao). More precisely, there exists a 
neighborhood U C X of u = such that U is separated into two disjoint 
open sets U x and by the stable manifold T\ of u = satisfying the 
following properties: 

(a) u = u? + u£+r x , 

(b) the transition in U x is jump, and 

(c) the transition in is continuous. The local transition structure is as 
shown in Figure \A.3[ 

(2) (A.l) bifurcates in U% to a unique singular point v x on X > Xq, which is 



an attractor such that for any <p € U% 



lim 

t — >oo 



\u(t, if) 



wher e u(t,tp) is the solution of (A.l). 



(3) (A.l) bifurcates on X < Xq to a unique saddle point v with the Morse index 
one. 

(4) The bifurcated singular point v x can be expressed as 

v x = -(^(AJ/ajVCfc-De! + o^ 1 '^). 

A. 3. Singular Separation. In this section, we study an important problem asso- 
ciated with the discontinuous transition of (A.l), which we call the singular sepa- 
ration. 



(1) An invariant set S of (A.l) is called a singular element 



Definition A. 2. 

if S is either a singular point or a periodic orbit. 
(2) Let Si C X be a singular element of (A.l) and U C X a neighborhood of 
Si. We say that (A.l) has a singular separation o/S at X = X\ if 

(a) (A.l) has no singular elements in U as X < Ai (or X > X\), and 
generates a singular element Si C U at X = X\, and 

(b) there are branches of singular elements S>, which are separated from 
Si for X > Xi (or X < Xi), i.e., 

lim max dist(x.Y.\) = 0. 

A special case of singular separation is the saddle-node bifurcation defined as 
follows. 



Definition A. 3. Let U\ G X be a singular point of (A.l) at X = X\ with u\ ^ 0. 
We say that (A.l) has a saddle-node bifurcation at (iii,Ai) if 

(1) the index of L\ + G at {u\, Ai) is zero, i.e., ind{—{L\ l + G),ui) = 0, 

(2) there are two branches T \(X) and T 2 (A) of singular points of (A.l), which 
are separated from U\ for X > Xi (or X < Xi), i.e., for any u\ <E r^(A) 
(i = 1,2) we have 

u\ —> Ui in X as X — > Ai, 

and 



(3) the indices of u\ £ Ti(A) are as follows 

' 1 



id{-{L x + G),u x ) = 



1 



ifu x e r 2 (A), 
tf« A eri(A). 
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Intuitively, the saddle- node bifurcation is schematically shown as in Figure |A.4| 
where the singular points in T^A) are saddle points and in r 2 (A) are nodes, and 
the singular separation of periodic orbits is as in shown Figure |A.5| 




Figure A. 4. Saddle-node bifurcation. 




Figure A. 5. Singular separation of periodic orbits. 



For the singular separation we can give a general principle as follows, which 
provides a basis for singular separation theory. 



Theorem A. 4. Let the conditions (A. 4-) and (A. 5) hold true. Then we have the 
following assertions. 



(1) If (A.l) bifurcates from (u, A) = (0, Ao) to a branch Y,\ of singular elements 



on A < A which is bounded in X x (— oo, Ao) then ( pl.i[ ) has a singular 
separation of singular elements at some (Eo,Ai) C X x (— oo, Ao). 
(2) If the bifurcated branch T.\ consists of singular points which has index —1, 
i.e., 



ind(-(Lx + G),u x ) 



-1 Vw A e E x , A < Ao, 



then the singular separation is a saddle-node bifurcation from some (u%, X\) G 
X x (— oo, A ). 



We consider the equation (A.l I defined on the Hubert spaces X = H, X\ = H\. 
Let L x = -A + \B. For L\ and G(-, A) : H 1 
symmetric, and 



H, we assume that A : Hi — * H is 



< Au,u >h> c \\ u \\'h 1 



/2' 



(A.7) 

(A.8) < Bu , U >H> C \\ U \\ 2 H 7 

(A.9) <Gu,u> H <- Cl \\u\\ p H 

where p > 2, c, C\, c 2 > are constants. 



c 2||"llff, 
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Theorem A. 5. Assume the conditions (A. 3), (A. 4) and (A.l)- \A.9jj , then (A.l) 
has a transition at (u, A) = (0, Aq), and the following assertions hold true: 

(1) If u = is an even-order nondegenerate singular point of L\ + G at A = Ao, 
then ( L4.i| ) has a singular separation of singular points at some (i(i,Ai) S 
H x (—00, A ). 

(2) If m = 1 ofirf G satisfies (A. 6) with a > if k=odd and a^O if k=even, 
then (A.l) has a saddle-node bifurcation at some singular point (ui,Ai) 
with Ai < Ao- 

A. 4. Transition and Singular Separation of Perturbed Systems. We con- 
sider the following perturbed equation of ( |A.1[ ): 

du 



(A.10) 



dt 



= (L X + S e x )u + G(u,\)+T"(u,\), 



X is a linear perturbed operator, 



where L\ and G\ are as in (A.l), S\ 

: X a — > X a C 1 nonlinear perturbed operator, and X a the fractional order 
space, < a < 1. Also assume that G\,T^ are C 3 on u, and 



(A.ll) 



< e, 



T £ (u,A) =o(||u|UJ. 



Let (A.4| and (A. 5) with m = 1 hold true, A) = G 2 (u, A) +o(||u|| j^), where 
G2O, A) is a bilinear operator, and 

(A.12) b=< G 2 (e,A ),e* >^ 0, 

where e£l and e* £ A* are the eigenvectors of L\ and corresponding to /3i(A) 
at A = Ao respectively. 

We now consider the transition associated with the saddle-node bifurcation of 
the perturbed system (A. 10 1. Let h(x,X) be the center manifold function of (A.l) 
near A 



(A.13) 



Ao- Assume that 

< G(xe + h(x, A ), A ), e* >= bix 3 + o(|cc| 3 ), 



where b\ 7^ 0, and e and e* are as in (A.12 1. 
Then we have the following theorems. 



Theorem A. 6. Let the conditions (A. 4) and (A. 5) with m = 1, and (A.13) hold 
true, and &i < 0. Then there is an e > such that if S x and T| satisfy (A.ll), 
then the transition of (A. 10) is either continuous or mixed. If the transition is 
continuous, then Assertions (2) and (3) of Theorem A. 2 are valid for ( A.10\ . If 
the transition is mixed, then the following assertions hold true: 

(1) ( A. 10) has a saddle-node bifurcation at some point (ui, Ai) G A"x(— 00, Aq), 
and there are exactly two branches 

Lf = {K A ,A)| Ai < \<\ £ Q + 6} i=l,2, 



separated from (u 1; Ai) as shown in Figure [A.4\ which satisfy that 

Vh\A)er^, Ai<A<A§ + <f, 
|«iHa- = vK,A)eL A . 



H\\x + 



lim 
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(2) There is a neighborhood U C X of u = 0, such that for each X with Ai < 
A < Aq + S and X 7^ Ag, U contains only two nontrivial singular points u A 
and n A of {A.10\ . 

(3) For each Ai < A < Aq + <5, U can be decomposed into two open sets U = 
t)£ + U£ with n t/ 2 A = 0, such that 

(a) ifXx < A < A§, 



S Ut 



u A € U^, 



u A e <9[/ A n<9[/ 2 A , 



wi£/i n = and n 2 6emj attractors which attract t/ A and \J 2 respec- 



tively, and 
(b) ifX e <X<X £ + 5, 

4 € 



ul e /7 2 \ 



o e ac/^ n <9[/ 2 \ 



wi£/i n A and u 2 feeing attractors which attract /7 A and £7^ respectively. 
(4) iVear (u, A) = (0, Ag),u A and u 2 can be expressed as 

u A = ai(A,e)e + o(|ai|), 



(A.14) 



a 2 {X,e)e + o(\a 2 \), 
0, 



lira ai (A, e) ■■ 
a 2 (A^,e) ^0, 



where e is as in {A. 13). 



Theorem A. 7. Assume the conditions (A. 4) and (A. 5) with m — 1, and (A. 13) 



with b\ > 0. The 



fh 



ere is an 



e > such that when and T| satisfy (A. 11), 



the transition of (A. 10) is either jump or mixed. If it is jump transition, then 
Assertions (1) and (3) of Theorem A. 2 are valid for (A. 10). If it is mixed, then 
the following assertions hold true: 

(1) (A. 10) has a saddle-node bifurcation at some point (mi,Ai) G ATx(Ag,+oo), 
and there are exactly two branches 

r A = {(u A ,A)| XI -5 <A< Ai} (1 = 1,2), 



5 < X < Ai} 
separated from (ui,X\), which satisfy 

o v(4,A)er; 



e < A < Ai, 



lira 



Kllx =0 



V(u A ,A) er A . 



(2) There is a neighborhood U C X ofu = Q, such that for each X with A g — 5 < 
X < Xi,U contains only two nontrivial singular points u\ and it A of (A. 10). 

(3) For every Ag — 5 < X < X\,U can be decomposed into three open sets 
U = U + Ui + U 2 with Ui tlUj = (ij^ j) such that 

(a) if Ag — 5 < X < Ag, then 

u = 0£ C/ \ v£ S <9t/ A n dU^(i =1,2), 

with u = oeing an attractor which attracts Uq and J7 A (i = 1,2) two 
saddle points with the Morse index one, and 

(b) if Ag < A < Ai, then 



n A e C/ A , 
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with Ui being an attractor which attracts and u% and u = being 
saddle points with the Morse index one. 
(4) Near (0, \®),Ui and can be expressed by 
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